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Abstract 

We perform both classical molecular dynamics simulations and beam model calculations to 
investigate the Young's modulus of kinked silicon nanowires (KSiNWs). The Young's modulus is 
found to be highly sensitive to the arm length of the kink and is essentially inversely proportional to 
the arm length. The mechanism underlying the size dependence is found to be the interplay between 
the kink angle potential and the arm length potential, where we obtain an analytic relationship 
between the Young's modulus and the arm length of the KSiNW. Our results provide insight into 
the application of this novel building block in nanomechanical devices. 
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I. INTRODUCTION 



In 2009, a new type of building block - kinked silicon nanowire (KSiNW) - was synthe- 
sized by Tian et al. in Lieber's group at Harvard University^. The growth direction of 
the KSiNW changes from < 112 > arm to < 110 >joint to < 112 > arm at the kink. Par- 
ticularly, the researchers are able to manipulate the arm length of the kink by controlling 
the growth time. Great experimental efforts have been devoted to investigating the growth 
mechanism of the kinks in silicon nanowires since thenar-. Besides these experimental works, 
Schwarz and Tersoff have proposed a theoretical model to interpret the growth mechanism 
of the KSiNW. In their model, the kink formation is due to the interplay of three basic 
processes: facet growth, droplet statics, and the introduction of new facets^. Stimulated 
by KSiNWs, several groups have synthesized kinks in other nanowires, such as 1^03 multi- 
kinked nanowires^, kinked germanium nanowires^, germaniumsilicon axial heterostructure 
with kinks^, and kinked ZnO nanowires^. While existing works for KSiNW mainly con- 
centrate on its growth mechanism, the study of elastic properties like Young's modulus is 
also important for its application in nanomechanical devices. 

The bulk silicon is an anisotropic mechanical material with the value of Young's modulus 
distributed roughly between 110 GPa and 180 GPa depending on the lattice direction.— ~— 
Hopcroft et.a/have shown that the Young's modulus in bulk silicon along < 110 > direction 
is about 45% higher than that in the < 100 > lattice direction.— In nano materials, due 
to large surface to volume ratio, the size effect^ and nanodefects^ have been found to be 
important for the mechanical properties such as the Young's modulus. For example, the 
AFM-measured Youngs modulus ranged from 93 to 250 GPa depending on the nanowire 
diameter.— 

In this paper, we apply both molecular dynamics (MD) simulations and finite element 
method to study the size- dependence of the Young's modulus in the KSiNW. We find that the 
Young's modulus is sensitive to its arm length; specifically, it decreases fast with increasing 
arm length. The underlying mechanism for this size- dependence is disclosed to be the 
competition between kink angle potential and the arm length potential. 
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II. LATTICE CONSTRAINT ON THE DIAMETER OF THE KSINW 



Fig. [TJ displays the configuration of a KSiNW. Panel (a) shows the growth axis of the 
KSiNW in a cubic lattice. Si, a 2 , and 0*3 are the three primitive vectors in the cubic lattice, 
with lattice constant a = 5.43 A. Lattice vectors i?i 2 = (2, 1, 1) and i? 34 = (1,2, —1) form 
two arms of the kink, while -R23 = (1, 1,0) forms the joint of the kink. Geometrically, the 
KSiNW is denoted by a pair of integers (n arm , nj oint ). The arm length is b = n arm i?i 2 and 
the joint length is rij oint i?23- 

From the schematic view shown in panel (b), the diameter of a KSiNW is restricted to 
some discrete values as determined by the lattice vector i? 23 , i-e d = 2r 2 3 = 2nj oint i? 2 3 = 
1.54njoint nm - This lattice constraint will have direct result on the configuration of the kink. 
In the experiment, the diameter d is pre-defined during the growth of the KSiNW-, which 
will likely deviate from the above discrete values. There are two possible consequences for 
the KSiNW to accommodate the pre-defined diameter. Firstly, the diameter of the joint 
part may become more difficult to be controlled in the experiment, if the diameter of the 
arm deviates from the discrete values. As a result, there will be a discontinuity at the kink 
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between < 112 > arm and < 110 >j i n t growth directions. r 
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the lattice constraint. The obvious discontinuity in Ref. [1 
Another possible way to break the lattice constraint is that the axial length of the joint 
becomes more difficult to be controlled, when the diameter of the arm deviates form the 
discrete values. As a result, the joint is not a perfect regular triangle anymore. Instead, 
the joint turns into a echelon-like structure, which is clearly shown in the Fig.2(c) of Ref. [H 
In KSiNW samples with diameter down to a few nanometers, the lattice constraint will 
play a more important role. Panel (c) shows four unit cells of a KSiNW with 1.54 nm in 
diameter and 2.66 nm in arm length, i.e (n arm , fijoint) — (2, 1). The unit cell is highlighted 
by a rectangular box. 

III. RESULTS AND DISCUSSION 

A. Results from MD simulation and finite element method 

In our MD simulation, the interaction between Si atoms is modeled by the widely used 
Stillinger- Weber empirical potential. The structure of the KSiNW is first relaxed without 



applying strain, yielding the minimum total energy Eq. After applying strain e to the 
system, the tensile structure is optimized with left and right ends fixed, which results in a 
total energy E{e). The axial Young's modulus is then calculated from the second derivation 
of the strain energy density {E(e) — E )/V with respective to the strain e in [0, 0.01]; V is 
the volume. It should be noted that the obtained Young's modulus is that of the overall 
kinked system (i.e. parallel to the unit cell in Fig. [TJ (c)), and not the Youngs modulus in 
the axial direction in each part of the nanowire. 

A sensitive size-dependence is observed for the Young's modulus as shown in Fig. [2] for 
three sets of KSiNWs with diameters 1.54, 3.07, and 4.61 nm. Only one unit cell is considered 
here, and we found that the Young's modulus does not depend on the number of unit cells 
in the KSiNW. The Young's modulus decreases with increasing arm length and can be fitted 
to an analytic function (solid lines) y = Ci/(feo + C2/&0) which can be well explained by the 
valence force field model as will be shown later. The two coefficients (01,02) are (216.5, 
0.23), (591.1, 1.45), and (881.7, 2.12) for these three sets of KSiNWs. It shows that the 
Young's modulus can be very small for KSiNW with very long arm length. This prediction 
can be readily verified experimentally, considering the arm length of KSiNW samples can 
be successfully controlled in the laboratory^. 

The inset of Fig. [2] shows the diameter dependence for the Young's modulus in KSiNWs 
with arm length 13.3 nm. The Young's modulus increases exponentially and saturates to a 
constant value of 151.0 GPa. Similar diameter dependence was observed in straight silicon 
nanowires.— Using the Stillinger- Weber potential, we have calculated the Young's modulus 
in a bulk silicon to be 159.0 GPa in < 110 > direction and 132.0 GPa in < 112 > direction. 
Both values fall in the experimental range of the Young's modulus in bulk silicon.— ~— The 
saturation value of 151.0 GPa in large diameter limit is sandwiched between the Young's 
modulus in bulk silicon along < 110 > and < 112 > directions. This result can be understood 
from the fact that the KSiNW is constructed by < 110 > (the arm) and < 112 > (the kink) 
silicon nanowires. 

We also calculate the Young's modulus of KSiNWs from the continuum theory based on 
the finite element method. We use the Euler-Bernoulli beam model22r— in our calculation. 
Finite element calculations are performed using the commercial ANSYS 12.0 package with 
2-node BEAM188 element. The only input parameters required by this package is the 
Young's modulus and Poisson ratio for bulk silicon from the above calculation based on the 
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Stillinger- Weber empirical potential. Fig. [3] shows the Von Mises stress distribution from the 
elastic beam model for two KSiNWs with arm length 13.30 nm and diameter 1.54 nm in (a), 
and arm length 26.60 nm and diameter 1.54 nm in (b). The applied tensile displacements 
are the same in both structures. The stress in the thinner KSiNW in (b) is roughly one 
quarter of that in the KSiNW in (a), because the Young's modulus and the strain in the 
thinner KSiNW in (b) are both about half of that in the KSiNW in (a). 

We note that the Euler-Bernoulli beam model neglects the shear strain, while this shear 
strain is considered in the Timoshenko beam model. We have also tested the Timoshenko 
beam model, and a similar size dependence of the Young's modulus is obtained while the 
value is overall smaller. An interesting phenomenon in Fig. [3] is the distinct stress concen- 
tration at kinks. The stress concentration may lead to the generation of a fracture around 
the kink in the KSiNW under large tension. This stress concentration phenomenon also 
affects the accuracy of the predicted Young's modulus. An important parameter in such 
calculation is the stress concentration ratio (SCT), i.e the ratio of the stress concentration 
length to the total length. The SCT is less than 30% in both KSiNWs shown in Fig. |3j This 
is important. It is because the Saint- Venant principle says that it is crucial for the SCT to 
be less than 30% for an accurate prediction of the Young's modulus.— 1 ^ According to the 
Saint- Venant principle,— a more accurate Young's modulus can be obtained from a smaller 
SCT, because the stress is distributed uniformly in more areas within the system in case of 
smaller SCT. From the comparison of these two KSiNWs in Fig. |3l the SCT decreases with 
increasing arm length. As a result, the Young's modulus from the finite element method 
should be more accurate for KSiNWs of longer arm length. Indeed, Fig. |2] shows a better 
agreement between the Young's modulus from the finite element method and the molecular 
mechanics approach for KSiNWs of longer arm length. 

Results for the set of KSiNWs of diameter 1.54 nm are shown by the dashed line in 
Fig. [2j The overall tendency from the finite element calculation agrees well with the MD 
simulation results. The Young's modulus decreases with increasing arm length. It should 
be noted that the finite element results always lie bellow the MD simulation results, which 
may probably due to the lose of part of the angle bending interactions in the beam model 
and also surface effects. As a result, the beam model actually mimics a system that is softer 
than the real material. The softness property of the beam model has been discussed in more 
detail in a recent work by Zhao et.al— The deviation between the finite element results 
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and the MD simulations is more pronounced for KSiNWs with smaller arm length, as the 
Euler- Bernoulli beam model is more accurate for a thinner rod structure with higher length 
to diameter ratio. There will be some systematic errors in the Euler- Bernoulli beam model 
if the simulated is very thick and/or short, since the higher order terms associated with the 
shearing deformation are ignored in this model. The surface induced diameter dependence 
of the Young's modulus from the finite element method is also shown in the inset of Fig. [2j 
Again, the general tendency from finite element method agrees with that from the MD 
simulation, and the elastic method also gives a lower value for the Young's modulus. 

We note that it is a usual technique to use the bulk Young's modulus of silicon as input 
for finite element modeling of mechanical properties of nanomaterials or nanostructures.— ^ 
Yet, there is some uncertainty in doing so, since the nanowires are anisotropic in atomic- 
scales. A different Young's modulus value may lead to an overall shift for the whole curve in 
the inset of Fig.[2j while the diameter dependence of the Young's modulus is kept unchanged. 
As a result, the size dependence of the Young's modulus from the finite element modeling 
still agrees with that from the molecular mechanics method, although the input Young's 
modulus for the finite element modeling varies. This size dependence is more important in 
the present work. In this sense, it is appropriate to use the bulk properties of silicon for the 
finite element modeling of the kinked silicon nanowires. 

B. Analysis based on valence force field model 

The above size-dependence for the axial Young's modulus can be understood in a valence 
force field model. Considering its zigzag configuration, the KSiNW can be simplified by 
a series of springs representing its arm length and kink angle as shown in Fig. HI The 
arm-spring has a force constant of kb, while the kink-spring has a force constant of ke- 
The potential of these two springs give the most important interaction for such a zigzag 
structure. Due to the translation invariance, we only need to consider the interaction within 
a single kink in the KSiNW. The kink angle bending (Vg) and the arm length stretching (V&) 
potentials are^ - — , 

Kot = V e + V b = ^- (cos 9- cos^o) 2 + § (6 - hf , (1) 
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where 9 is the kink angle and b is the arm length. Variables with subscript correspond to 
their values at equilibrium configuration without strain. For uniaxial strain applied in the 
horizontal direction, Fig. H] shows a relationship among geometrical variables 
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b sin 



(2) 



There is only one independent degree of freedom in the tensile KSiNW. We choose the kink 
angle 9 as the free variable. From the energy minimum condition dE/dO = 0, we have 



kf> (cos 9 — cos 9 n ) sin 9 sin 3 - 



+ ^hc ( c - & o sin - J cos - = 0. 



(3) 



The solution of Eq. ([3]) yields the equilibrium structure of the KSiNW under mechanical 
strain e. For small strain, the kink angle is also small, so it only slightly deviates from its 
equilibrium value, i.e 9 = 9 + 5g. Applying standard perturbation theory to Eq. ([3]) up to 

(8 e ) 2 , we get 



a6$ + /35 e + 7 = 0, 



(4) 



where the three coefficients are: 

a 

P 



hc(fb -± 
he (Jfto - 4 c) - ^k g 



(5) 



1 = k b c(\c-^b ). 



The solution of the equilibrium equation is then obtained analytically, 



Si 



2a 



(6) 



It can be shown that the elastic properties are dominated by the first-order term in the 
Taylor expansion of the angle variation Sg in terms of e. We thus expand the three coefficients 
in Eq. (jSJ) to the Taylor series of e: 



a ~ «o + cx\e, 
P ~ Po + Pie, 
1 ~ 7o + 7ie, 



(7) 
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where 



a 

7o : 



Ikf,2 
64 i^b u o, 



V3t, 13 9V3L, 

0. 



(8) 



Inserting these equilibrium coefficients into Eq. (J6]), we obtain the angle variation without 



strain: 8a 



\M 

a 



and 8g = 0. Obviously, only 8 e is the physical solution, because the angle 



variation should vanish without strain. For the first-order term, we have 



The angle variation is: 



[ «1 = 


dot! h.2 

de l £=0 ~~ 32 ft6W 0> 




Pl = 


de 1^=0 — 4 ft feO , 


(9) 


k 7i = 






8 e = 


7i 1 


(10) 


\/3 \ C f + ^rj^ 



where rj = According to this formula, the kink angle variation is determined by the 
competition between kink angle and arm length potentials. It is smaller for KSiNW with 
larger kg, or smaller k b , which is physically correct. A rigid arm (kb — > +oo) results in the 
largest angle variation of 8g = 2\^3e. 

Inserting the kink angle variation into the total strain energy in Eq. (fTTh and do linear 
approximation, we obtain the elastic strain energy density, 



Ah 



'0 



^kf>8p 
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where A is the cross section of the KSiNW. The Young's modulus is 
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where coefficients Ci=a and c 2 = 

Eq. (TTTT) is exactly the fitting function for the MD simulation results in Fig. |2j This 
coincide verifies the spring representation of the KSiNW. We note that the deviation between 
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MD simulations and the valence force field model in the figure for longer arm length is 
probably due to nonlinear effects resulting from larger kink angle variation for longer arm 
length, which is omitted in our derivation. From the above analytic derivation, it is clear 
that the Young's modulus is determined by the competition between the kink angle bending 
potential and the arm length stretching potential, corresponding to the two terms in the 
denominator of Eq. (ITT]) . For KSiNW with bo < ^Joi, the Young's modulus increases with 
increasing arm length; while the opposite phenomenon should be observed in KSiNW with 
fro < v/c2- Our MD simulation results in Fig. [2] show that these KSiNW in our study fall in 
the latter case. The size- dependence of the Young's modulus is governed by the coefficient 
C2, which is determined by the ratio of the two force constants kg and kf,. Hence, our analytic 
formula provides valuable information for tuning the behavior of the Young's modulus by 
using various materials with different kg or kf,. It is interesting that C2 is independent of 
the diameter of the KSiNW, so the decreasing behavior of the Young's modulus is diameter 
independent, which is indeed observed in the MD simulation results. 

Using the fitting parameters (01,02) from our MD simulations, we can extract the two 
force constants {k b , kg) to be (103.8, 265.3), (179.8, 2897.3), and (412.8, 9723.9) for KSiNWs 
with diameters 1.54, 3.07, and 4.61 nm, respectively. Dimensions of these two force constants 
are [kb}=eV/A 2 and [fcg]=eV. Both force constants fc& and kg are independent of the arm 
length, but increase rapidly with increasing diameter. These two force constants are effective 
values for the whole KSiNW system, so it can be quite different from the force constants 
in the potential for an actual chemical bond. For instance, kg for KSiNWs with diameter 
1.54 nm is about 20 times larger than that describing C-C-C angles in carbon nanotubes^, 
while kb is half of that in the C-C bonds in carbon nanotubes^i^. From the energy point of 
view, each silicon atom in the system is connected to its neighbors by small springs within 
the linear approximation, kb and kg are actually the summation over these small springs 
between silicon atom pairs. For KSiNW with larger diameter, there are more small springs 
that are in parallel connection, so the two effective force constants k b and kg are larger. For 
KSiNW with longer arm length, there is an increasing number of small springs that are in 
series connection, which shouldn't affect the effective force constants. As a result, kb and kg 
is sensitive to the diameter of the KSiNW, while is independent of its arm length. 
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IV. CONCLUSION 



To summarize, we first point out a lattice constraint on the diameter of the KSiNW. 
Then we perform extensive MD simulations and finite element method to study the Young's 
modulus of the KSiNW. The Young's modulus is found to decrease rapidly with increasing 
arm length. The size sensitivity is explained by a valence force field mode, where we derived 
an analytic formula for the axial Young's modulus of the KSiNW with different arm length, 
i.e Y = c\/{b + C2^). We show that the size- dependence of the Young's modulus in KSiNW 
is governed by two different types of potentials stored inside its zigzag configuration, which 
lead to the two terms in the denominator. 
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(c) 




FIG. 1: (Color online) The configuration of a KSiNW. (a) shows the growth axis from < 112 > arm 
to < 110 >joint to < 112 >arm- (b). A schematic geometry, disclosing a lattice constraint on 
the diameter d, i.e d/2 = r23- (c). Configuration of a real KSiNW with structure parameters 
("arm) ^joint) = (2, !)■ The rectangular box (blue online) highlights a unit cell of the KSiNW. 
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FIG. 2: (Color online) The size effect on the axial Young's modulus. The MD simulation results 
for the Young's modulus versus arm length in KSiNW with diameters 1.54 nm (black points), 
3.07 nm (blue triangulars), and 4.61 nm (red squares). Simulation results are compared with the 
valence force field model (solid lines) and the elastic beam model (dashed lines). Inset: the Young's 
modulus versus diameter for KSiNW with arm length 13.30 nm, where calculation results (points) 
are fitted to an exponential function (solid line) and compared with the elastic beam model (dashed 
line). 
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FIG. 3: (Color online) Von Mises stress (in unit of 10 10 Pa) distribution in KSiNWs from the elastic 
beam model calculation. The applied tensile displacement is the same in (a) for the KSiNW with 
arm length 13.30 nm and diameter 1.54 nm, and (b) for the KSiNW with arm length 26.60 nm 
and diameter 1.54 nm. Note the stress concentration at kinks. Attention that the longer nanowire 
in (b) has been scaled by an overall factor of 0.5. 




FIG. 4: (Color online) The KSiNW is simplified by a valence force field model with spring constant 
kf, and angle bending constant kg. A single kink (inside the box) is highlighted for its equilibrium 
structure without strain (solid, black online) and under mechanical strain (dotted, blue online). 
Geometrical variables are related to each other by c = b sin | . 
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